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Abstract 

We study the perturbative behaviour of topological black holes. We calculate both 
analytically and numerically the quasi-normal modes of scalar perturbations. In the 
case of small black holes we find discontinuities of the quasi-normal modes spectrum 
at the critical temperature and we argue that this is evidence of a second-order phase 
transition. 
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1 Introduction 

The knowledge of the spectrum of the quasi-normal modes (QNMs) is a powerful tool in 
the study of the late time behaviour of black holes. The QNMs are the complex frequencies 
by which a black hole responds if it is initially perturbed, and they do not depend on the 
details of the initial perturbation but rather on the intrinsic features of the black hole itself. 
The radiation associated with these modes is expected to be seen with gravitational wave 
detectors in the coming years, giving valuable information on the properties of black holes. 

The QNMs of black holes in asymptotically flat spacetimes have been extensively stud- 
ied and their spectrum was computed numerically and in many cases also analytically 
(for reviews, see |U [2j). The advances in string theory and mainly the Anti-de Sitter - 
conformal field theory (AdS/CFT) correspondence has renewed the interest of computing 
the QNMs of black holes in asymptotically AdS spacetimes [21 HJ El El El E] • Recent re- 
sults from the Relativistic Heavy Ion Collider [H] show that a thermal quark-gluon plasma 
(QGP) is formed which is strongly coupled. The AdS/CFT correspondence provides the 
connection between the QGP and string black holes. In [10] the quasi-normal modes of 
AdSs-Schwarzschild black hole have been calculated and it was shown that they provide 
a dual description of the fluctuations of the QGP. Recently, in a more phenomenological 
approach, the AdS/CFT correspondence was applied to condensed matter physics [TTJ. It 
was shown [12] that fluctuations of the metric and the background electromagnetic field 
determines the conductivity of the boundary theory. 

According to the AdS / CFT correspondence, a large static black hole in AdS corresponds 
to an (approximately) equilibrium thermal state in the CFT [TJ]. In ref. [I] it was shown 
that the QNMs for the scalar perturbations of large Schwarzschild-AdS black holes scaled 
with the temperature and it was argued that the perturbed system in the dual description 
will approach to thermal equilibrium of the boundary conformal field theory. However, 
when the black hole size is comparable to the AdS length scale there is a clear departure 
from this behaviour. It was then conjectured that this behaviour may be connected with a 
Hawking-Page phase transition [TBI [TB] which occurs when the temperature lowers. 

For small black holes the behaviour of QNMs is not very well understood. For the case 
of Schwarzschild-AdS black holes the QNMs do not scale linearly with the temperature any 
more and their connection with the boundary conformal field theory is not clear. In [T71 [T8] 
we studied the QNMs of electromagnetic and gravitational perturbations of topological 
black holes (TBHs) in AdS spacetime coupled to a scalar or an electromagnetic field. For 
small black holes, compared to the length scale of the AdS space, we found that the QNMs 
behave quite differently from the QNMs of large black holes. Near the critical temperature 
purely dissipative modes appeared in the spectrum and the QNMs change slope around 
that point. We attributed this behaviour to a second order phase transition of the charged 
TBH towards the AdS vacuum solution. 

In the literature there is a discussion of possible connections between the classical and 
thermodynamical properties of black holes [19]. In particular the question whether the 
knowledge of the QNM spectrum can give information about thermodynamical phase tran- 
sitions in a wider class of black holes has recently gained considerable interest. It was 
suggested in [20] that the Dirac and Rarita-Schwinger perturbations are related to ther- 
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modynamic phase transitions of charged black holes. This has been criticised in [21], the 
argument being that the relation between the QNMs and the phase transition had not 
been properly formulated. In this direction, the discontinuities observed in the heat ca- 
pacity of charged Kaluza-Klein black holes with squashed horizons, were connected with 
the quasi- normal spectrum [22]. Further evidence of a non-trivial relation between the 
thermodynamical and dynamical properties of black holes was provided in [23J. 

In this work we calculate both analytically and numerically the QNMs of scalar pertur- 
bations of topological- AdS black holes in d — 4, 5 and 6 dimensions. In all dimensionalities 
considered, we find a discontinuity of the QNMs spectrum at the critical point. This pro- 
vides further evidence of a second order phase transition at the critical point observed 
in [TTl [18] . In our numerical calculations we used the method developed in [I]. We conjec- 
ture that this method is problematic for large n because it breaks down, and a regularization 
scheme is proposed. 

Our work is organized as follows: in section 2 the analytical calculation is presented, 
section 3 contains comments on the numerical method used, and in section 4 the numerical 
results may be found. Finally the conclusions are presented in section 5. 



2 Analytical Calculation 

Analytical results for the QNMs of topological black holes are only known for the 
massless case. In four dimensions they were calculated in [30J while a generalization to 
d-dimensions have been presented in [31]. For the general case early numerical results have 
been discussed in [3] and in [29]. To study the phase transition in topological black holes 
the QNM spectrum has to be known, at least near the critical point. Therefore in this 
section we calculate analytically the QNMs of scalar perturbations of topological black 
holes in various dimensions. At the critical point the mass goes to zero so we reproduce the 
known exact results. Also we discuss the asymptotic form of QNMs for large black holes. 
In five dimensions we present analytical results expressible through the Heun function. The 
analytical results will give us a clear indication of a second order phase transition around 
the critical point which will be substantiated by numerical results presented in the next 
section. 



We consider the bulk action 

1 



16ttG 



(rf-l)(d-2) 
R+ f2 



(2.1) 



in asymptotically AdS^ where G is the Newton's constant and / is the AdS radius. The 
presence of a negative cosmological constant (A = — < ^—^ z ^-) allows the existence of black 
holes with a topology IR x E, where S is a (d — 2)-dimensional manifold of constant negative 
curvature. These black holes are known as topological black holes. The simplest solution 
of this kind reads 

ds* = -f(r)dt* + rW , f(r) = r 2 - 1 - ^- , (2.2) 



f(r) ' r 
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where we employed units in which the AdS radius is I = 1 and da is the line element of E. 
The latter is locally isomorphic to the hyperbolic manifold H d ~ 2 and of the form 

E = H d - 2 /T , rcO(d-2,l), (2.3) 

where T is a freely acting discrete subgroup (i.e., without fixed points) of isometries. The 
geometry of the topological black holes as well their basic properties have been studied 
extensively in the literature |24j-[33j. 

The configurations (12.21) are asymptotically locally AdS spacetimes. It has been shown 
in [31] that the massless configurations where S has negative constant curvature are stable 
under gravitational perturbations. More recently the stability of the TBHs was discussed 
in [SB] . 

We shall calculate analytically the quasi-normal modes of scalar perturbations of these 
black holes following the method discussed in [IB] . A numerical approach will be discussed 
in section [31 

Massless scalar perturbations obey the radial wave equation 

(r*- 2 f(r)&)' + - \ 2 J $ = . (2.4) 



By defining 



d-2 



$ = r -— m (2.5) 

the wave equation can be cast into a Schrodinger-like form, 

--fY + V[r{n)}V = 2 ^ , (2.6) 
in terms of the tortoise coordinate defined by 

dn 1 (2.7) 



dr f(r) 
and the potential is given by 



d(d-2) - f (d-2) 2 fi\ , 2 _ C 2 fd-3\ 2 



V{r) = f{r){^^ + ^ + ^^\ , k 2 = e+r~^) • (2.1 



For later convenience we define A = k 2 — 3/4. 

To obtain analytic expressions for the quasi-normal frequencies, it is convenient to 
introduce the coordinate 

u=( r -±Y. (2.9) 



r 

The wave equation (12.61) becomes 



Au z ' 2 f(u)(u z/2 f(u)m')' + [Cj 2 - V)V = , cu = — , (2.10) 

r + 
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where prime denotes differentiation with respect to u and we have defined 



+ 



and 



1 f 2\i ^3 
— — — — u 2 

1 H— 1 ' 

u r\ r\ 

f(u) t — - + Am + -r^u— 



2r\ 



d-i 



A 



' 4- 



(2.11) 



(2.12) 



2.1 Exact solution in d = 5 

In five dimensions (d = 5) we may solve the wave equation in terms of a Heun function and 
use the latter to determine the spectrum exactly albeit numerically. By factoring f(u), 

/(„) = (!-„) (1 + 1-^ 

we may obtain an exact solution to the wave equation in terms of a Heun function, 

i V 1 ~ 1/r 



*( u ) = u~ 3/4 ( u + 



1 



2(2-l/r 2 ) 



(1-u) 2(2 - 1/r + ) Heun(a, g, a, /3, 7, <?,!-«) (2.13) 



1 - l/r\ 

The Heun function obeys the equation 

z(z-l)(z-a)F"-{(-a-(3-l)z 2 + ((5+-f)a-5+a+f3+l)z-a-f}F'-l-af3z+q]F = (2.14) 
and the various constants are 



1 



1 - 



1 ' 



Q 



Cj\2 



2(1 -jr){2-jrj 
3 



^r(^ 2 + 1 )(2-i) 2 + ^(2--i) 



+ 



2r\ 



)+^-2 + i)(l-i) 3/2 



a = (3 



i+ A /l-4-)a) 



2(2-4) 2- 
It behaves nicely at the horizon (m — > 1). Requiring \l/(0) = yields the constraint 

Heun(a, q, a, (3, 7, 5, 1) = , (2-15) 

which may be solved for uj to obtain the quasi-normal frequencies of scalar modes in five 
dimensions. 
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2.2 Large Black Holes 

For large black holes, we let a), A — > 0, r + — > oo. The wave equation becomes 

M 3 / 2 [(u^ - „</* W - (d ~ 2)[d+ 1 ( 6 d ~ 2) ^U = . (2.16) 
Two linearly independent solutions are 

d-2 d-2 / d-l\ 

m = u~— , In ( 1 - tt~ ) . (2.17) 



Both are unacceptable as they diverge at u = 0,1, respectively. It follows that the only 
possible frequencies are those which diverge as r + — > oo. 

In five dimensions, the lowest mode may be deduced from the exact constraint in the 
limit r + — > oo keeping £ and Co fixed (and therefore A — > 0, uo ~ r + ). In this limit, 

= [2(« -!) + («- 2)g]g = = _<£+!£ _S , = _ x _ 

' y 8 P 4 ' ' 2 ' 

The lowest frequency in five dimensions is found to be 

Cj = 3.12-2.75? (2.18) 

in agreement with earlier numerical results [4J. Therefore, for a large black hole, u ~ r + 
and QNMs do not contribute to the macroscopic (hydrodynamic) behaviour of the gauge 
theory fluid on the AdS boundary. 

2.3 Near the Critical Point 

At the critical point (r + = 1, \i = 0), the wave equation reduces to 

J 1 d{d-2) 



1-u Au 
The solution which is well-behaved at the boundary is 



-A 



* = . (2.19) 



»(„) = (!- u) ^ F( ^J!i, - + !,„). (2.20) 

At the horizon, w — > 1, it behaves as 

*(«) ~ A(l - + - u)^ 2 , At = d+i r( f y )r( ^| f± . • (2.21) 



Demanding *4._ = 0, we deduce the quasi- normal frequencies 



u n = ±£-i [ 2n + ^J: ) , n = 0,l,2,... (2.22) 
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which have finite real part (except in the special case £ = 0). This result agrees with the 
quasi-normal frequencies of massless topological black holes which they were calculated 
analytically for the first time in [SO] . 

The hypergeometric functions corresponding to the eigenf unctions are polynomials. Ex- 
plicitly, 



*o(«) 



A (i- M )-w<mvv/4 

A n _ u )-«/2-(d+B)/V/4 



d+1 



(2.23) 



etc. They are orthogonal under the inner product (no complex conjugation!) 

" 1 du 



(n\m) 



o mV2(i _ w ) 



(2.24) 



defined by appropriate analytic continuation of the parameter £. To normalize them 
((n|n) = 1), choose 



42 



r(-«0 



(i+lN-p/d+lN ' 



-)r( 



(d + i)r(-i-«o 



(2.25) 



etc. 



Moving away from the critical point, the frequencies shift by 

1 (n\H'\n) 



5uj r 



2u) n (n\n) 



where 



(2.26) 



(2.27) 



and we applied standard first-order perturbation theory. 
We obtain 



5uj 



-i (r+ - 1) 



—i (r_i_ — 1 



2 d " 2 r(f )r(-^ - l)[-2if(d + 1) - d 2 + Qd - 17] 



(2.28) 



For small £, there is no change in the imaginary part at first order. Below the critical point 
(r + < 1), the change in the real part is negative (5u n < 0) and the real part decreases. 
There are critical values of £ (determined by £ + 8uj n ~ 0), 



2 d r(f) 



v^r( 



l-r+, 6 



, 2 ci - 2 r(f)[d 2 -6rf+17] 



v^r( 



(2.29) 



below which the corresponding mode does not propagate (purely dissipative mode). It 
turns out that for £ < £ n pairs of purely dissipative modes emerge. 
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Above the critical point (r + > 1), 5u n > and the real parts of the modes increase. 
The modes do not become purely dissipative for any value of £. 

Also notice that below the critical point, 8u n increases with n, therefore the real part 
decreases with n (positive slope) whereas above the critical point we obtain a negative slope 
for propagating modes. 

These results will be discussed further along with a comparison with numerical results 
in section HI 



3 Numerical Calculation 

In this section we present the numerical calculation of the QNMs of scalar perturba- 
tions of small topological black holes following the method developed by Horowitz and 
Hubeny jl]. 

After performing the transformation \l/(r) = ^e" 1Ijr *, the wave equation becomes 



d 2 ^(r) ( df(r) \ d^ w (r) 



dr 2 \ dr J dr 

The change of variables r = 1/x yields an equation of the form 



s(x) 



(x — *£+) 



dx 2 



+ t(x) 



dx 



+ u(x)i/j uj (x) = 



where x + = l/r + and s(x),t(x) and u(x) are given by 



k 

t( x ) = ^2t k (x-x+) k , 

k 

u ( x ) = ^2u k (x - x + ) k . 

k 

Expanding the wave function around the (inverse) horizon x + , 

oo 

i>w{x) = ^ a «M( x ~~ x +) n ' 



(3.1; 



(3.2) 



we arrive at a recurrence formula for the coefficients, 



. n— 1 

a n (uj) = -— — y a m (uj) 

' m—n— 4 



m,n—m 



m(m - l)s„_ m + mt„_ m + u r , 



(3.3) 
(3.4) 
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We note that the few coefficients a m (u) with negative index m which will appear for n < 2 
should be set to zero, while a,o(u) is set to one. Since the wave function should vanish at 
infinity (r — ► oo, x = 0), we deduce 

oo 

^(0)=^a n (o;)(-x + )" = 0. (3.5) 
o 

The solutions of this equation are precisely the quasi-normal frequencies. 
The change of variables 

Ct n = (X n X _j_ , Pm.n—m = Prn,n—m%-\- 771(771 l)s n _ m -|- Ttvt n _ m -\- 1t n _ m (3.6) 

transforms eqs. (13.31) and (13.51) into 

^ 71—1 

^.0 m =n-4 

and 

oo 

^ w (0) = 53a n (£ 1 ;)(-l) n = 0, (3.8) 
o 

respectively. 

To get a feeling about the solutions of the recurrence relations we consider the limit 
of very large n. In this limit we may keep just the terms of order n 2 , so relation (13. 6p is 
approximated by 

Pm,n—m — 171 S n — m (3.9) 

and the recurrence relation reads 

^ 71 — 1 

S ° m=n-N 

Notice that it has constant coefficients, so it can be solved analytically. For simplicity let 
us restrict attention to the case N = 4. The solution of the recurrence relation reads 

f^)Tl f^tTh f^Th 

It should be emphasized that the coefficients K m (m = 1,2,3,4) do not depend on n. 
Explicitly, 



rn 



where 



Am = s * (A 3 *r 6 m + A 2 *r m + A 1 * r m + A ) + s t * (A 2 * r m + A x * r m + A * r m ) 

+s 2 * (A 1 *r 3 m + A * r 2 J + s 3 * (A * r 3 J , 
B m = -(4 * s 4 * r m + 3 * s 3 * r m + 2 * s 2 * r m + si) . 
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r m are the four roots of the equation 

s + S\r m + s 2 r 2 m + s 3 r^ + s 4 r^ = , (3.10) 

respectively and Aj., k = 0, 1, 2, 3 are the values of the first four coefficients which are 
necessary to completely specify the solution. 

For large n, the sum is dominated by the root which has the smallest absolute value. 
If all roots have an absolute value which is bigger than one, the coefficients will go to zero 
for large n and the method is expected to work. However, if the absolute value of one of 
the roots is smaller than one, we run into problems. In our case, one of the roots is equal 
to -1. 

To see this, we report expressions for various coefficients in table [TJ Using these, one 



i 


Si 


u 


Ui 





x + (x\-3) 


—3x + + x+ + Hujx\ 





1 


x + (Ax\ - 9) 


-12x + + 6x% + Aiujx 2 + 


-3x + +d-a4 


2 


x+(6x% - 10) 


-18x+ + 12x 3 + + 2iujx 2 + 


-3a:+ + (| - 2e)x\ 


3 


x + (Ax\ — 5) 


-12x+ + 10x\ 


-3x + + (±-e)xi 


4 


x+(x^_ — 1) 


— Sx-^ -f- 3x^_ 


—x + + x+ 



Table 1: Coefficients used in the numerical method (eq. (13. 6p ) 
may show that 

Pn,0 + Pn-2,2 + -Pn-4,4 = Pn-1,1 + Pn-3,3 ■ 

Notice that this equality holds exactly in d = 4. For d = 5, 6 the two sums differ by an 
amount which is of order 0(1) so it is negligible for large n. In particular one may show 
the relation 

So + ^2 + S4 = Si + S3 , 

which implies that —1 is a root of eq. (I3.10p . as advertised. 

Thus our case is inconclusive and one should examine - corrections to settle the issue 

n 

of convergence. We expect that | | = l + -.lf/i<0 the sequence will converge to zero 
and the series can be finite. Otherwise there is a problem with convergence. 

To ensure convergence and obtain numerical solutions of the equations 3?[^(0)] = 
0, ^[^(0)] = , we shall introduce a regularization scheme as follows: rather than 
imposing the condition that the wave function vanish at r = 00, we impose the condition 
that it vanish at r = R = « 1. This translates to x = X = ex + . Then eq. (13.51) is 
replaced by 

00 00 
VUX) = J>nH(* - ^+) n = $> n M(-ir(l - e) n = . (3.11) 



Asymptotically, a n (uj) = (— l) n A(uj) which solves the recursion relation and is confirmed 
by numerical results. Without the reqularization parameter (e = 0), this leads to problems 
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Figure 1: Scalar QNMs in four dimensions at r + = 0.95, £ = 5.0. 

because the left-hand side of (13. 81) diverges. With a finite e, the series is asymptotically 
geometric and converges. 

If e is small enough, the error introduced is negligible. However, choosing too small a 
value for e is not calculationally efficient as it may call for too many terms in the series (so 
that the individual terms become small enough). We have used e = 0.01 for d = 4, d = 5, 
and e = 0.02 for d = 6. With these values our results agree with the analytical results for 
r + ps 1. 

4 Results 
4.1 d = 4 

We start with the case r + = 0.95 (below the critical point), £ = 5.0, depicted in figured) 
We observe that there exists a finite number of QNMs and positive slope. Similar results 
may be seen in figure [2j corresponding to r + = 0.99 and £ = 5.0. In the latter case the real 
parts do not change as much between successive modes. Going to smaller values of £ the 
most important changes are that the region of propagating modes becomes smaller, while 
the real parts take smaller values. This is apparent upon inspection of figure [3] (for £ = 3.0) 
and figure H] (for £ = 1.0). In the latter case purely dissipative modes appear. There exist 
no modes at all for imaginary parts absolutely larger than 35.. The purely dissipative 
modes disappear for larger values of £; for £ = 2 it is difficult to detect them. This fact 
is consistent with the analytic result that there exist critical values of £, above which only 
propagating modes exist. Thus, for small enough values of £ we have several propagating 
modes with relatively small absolute values of their imaginary parts and positive slope, 
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Figure 2: Scalar QNMs in four dimensions at r + = 0.99, £ = 5.0. 
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Figure 3: Scalar QNMs in four dimensions at r + = 0.99, £ = 3.0. 

followed by a finite number of purely dissipative modes. For larger values of £ only a finite 
number of propagating modes exists. 

To illustrate the features of the case with r + > 1.0, (above the critical point), we give in 
figure [5] the results for r + = 1.10, £ = 1.0. The slope is negative right from the beginning 
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Figure 4: Scalar QNMs in four dimensions at r + = 0.99, £ = 1.0. 



and no purely dissipative modes emerge. 
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Figure 5: Scalar QNMs in four dimensions at r + = 1.10, £ = 1.0 

Another set of results refers to the lowest modes, in particular the behaviour of their 
real parts as functions of £. The analytical results derived from eq. (12.281) along with the 
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numerical results are shown in figure [6j The agreement is good and it improves for large 
values of £. In particular, values of £ where these two real parts vanish compare quite well 
with the analytically computed values. 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 6: The analytical results (eq. (j2.28|) ) for the real parts of the two lowest QNMs 
uj and u)\ vs £ for d — 4 and r + = 0.999 are represented by lines, while their numerical 
counterparts by points. The upper (lower) curve and the corresponding points represent 
(^i)- 



In figure [7] we show the dependence of the critical values £o an d £i on the horizon r + . 
This figure depicts the analytical results expressed by eq. (12.291) . as well as the corresponding 
numerical results. The agreement between the two approaches is quite good. 



0.6 



Figure 7: The points represent numerical estimates of the critical values £o and £i as 
functions of the horizon for d = 4, while lines give the corresponding analytical estimates 
(eq. (12.291) ). The larger (lower) values represent £ (£i). 



Agreement between analytical (eqs. (I2.22j) and (12.281) ) and numerical results for the two 
lowest QNMs is further demonstrated in table [2j Let us comment that the two biggest 
critical values for £ in this case occur for r + = 0.95 and read £o — 0.50, £i = 0.76, well 
below the values of £ considered here. The agreement is better as we approach the value 
r + = 1.0, as is evident in the last four rows of table EJ 
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, .anal 


, ,num 


, .anal 


, ,num 


1.0 


0.95 


0.64 - 2.23i 


0.80 - 2.36i 


0.11 - 3.89i 


0.42-4.17i 


3.0 


0.95 


2.80 - 2.32i 


2.90 - 2.42i 


2.51-4.08i 


2.75 - 4.30i 


5.0 


0.95 


4.85 - 2.36i 


4.93 - 2A4i 


4.63-4.16z 


4.81 - 4.34i 


1.0 


0.99 


0.93 - 2.44i 


0.96 - 2A7i 


0.82 -4.38i 


0.91 - 4.44i 


5.0 


0.99 


4.97 - 2A7i 


4.99 - 2.49z 


4.93 - 4.43* 


4.96 - AA7i 


5.0 


0.999 


4.9971 - 2.4972i 


4.9985 - 2.4986i 


4.9926 - 4.4932i 


4.9963 - 4.4965z 


5.0 


0.9995 


4.9985 - 2.4986i 


4.9993 - 2.4993i 


4.9963 - 4.4966i 


4.9981 - 4.4982z 


5.0 


0.9999 


4.9997 - 2.4997i 


4.9998 - 2.4998i 


4.9993 - 4.4993i 


4.9996 - 4.4996z 



Table 2: Comparison of analytic vs numerical results for the two lowest QNMs at d = 4 



Figure 8: Imaginary part of the lowest purely dissipative QNM in four dimensions vs r + 
for £ = 1.0 (crosses) and £ = 2.0. 

Finally we depict in figure M the behaviour of the lowest purely dissipative QNM as 
a function of r + for £ = 1.0 and £ = 2.0. It tends to — oo as the horizon approaches the 
critical value r . = 1. 



4.2 d = 5 and d = 6 

Only quantitative changes appear in the five- dimensional case. We depict sample results 
in figure M Comparing against figure El let us only note the quantities |3?(u;)| change 
more between successive modes. We also note that in five dimensions the propagating 
modes coexist with purely dissipative modes in the same range of imaginary parts. Of 
course, there are no propagating modes with sufficiently large values of |Sy(a;)|, i.e., the 
propagating modes disappear and only dissipative modes survive. 

As one moves to d = 6, numerical difficulties become more severe and QNMs with too 
large imaginary parts become increasingly harder to explore. We show sample results in 
figure [TOj It seems that we encounter a qualitative change. The propagating QNMs have 
positive slope in the region of (absolutely) small imaginary parts, then there is a number 
of purely dissipative modes, while for even larger imaginary parts a new branch appears, 
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Figure 9: Scalar QNMs in five dimensions at r + = 0.99, £ = 3.0 
which has negative slope. 
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Figure 10: Scalar QNMs in six dimensions at r + = 0.9999, £ = 1.0. 

We also study the real parts of the lowest modes for d = 5 as functions of £ and present 
both the analytical results derived from eq. (I2.28j) and their numerical counterparts in 
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figure [ITJ The numerical results lie a bit above the prediction and the agreement is not so 
good as in the d = 4 case. This feature is even more pronounced for d = 6. The results 
are depicted in figure [T2"l It is evident that the most sizeable discrepancies occur for uj x . 
In particular, the numerical results appear more steep and the corresponding critical value 
for £i is substantially bigger than the analytical prediction. These results show that first- 
order perturbation theory is inadequate in this case - one needs to include higher-order 
perturbative effects. 
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Figure 11: The analytical results (eq. (I2.28P ) for the real parts of the two lowest QNMs 
u and u\ vs £ for d = 5 and r + = 0.999 are represented by lines, while their numerical 
counterparts by points. The upper (lower) curve and the corresponding points represent 
(^i). 
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Figure 12: The analytical results (eq. ( 12.28B ) for the real parts of the two lowest QNMs 
ujq and u% vs £ for d = 6 and r + = 0.999 are represented by lines, while their numerical 
counterparts by points. The upper (lower) curve and the corresponding points represent 

(^l)- 

In figure [131 we show the dependence of the critical values £o and ^ on the horizon r + 
in five dimensions. This figure depicts both the analytical results expressed by equation 
( I2.29P and the corresponding numerical results. Similar results are presented in figure [TJ] 
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Figure 13: The points represent numerical estimates of the critical values £o & n d £1 as 
functions of the horizon for d = 5, while lines give the corresponding analytical estimates 
(eq. (12.291) ). The larger (lower) values represent £o (£i)- 
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Figure 14: The points represent numerical estimates of the critical values £o and (i as 
functions of the horizon for d = 6, while lines give the corresponding analytical estimates 
(eq. (12.291) ). The larger (lower) values represent £o (£i)- 

for six dimensions. The agreement between analytical and numerical results is satisfactory 
in all cases. 

Table [3] contains the comparison of analytical versus numerical results for the two lowest 
QNMs at d — 5. The agreement is fairly good. Let us comment that the two critical values 
for £ read in this case: £o = 1.09, £i = 1.90 for r + = 0.95 and and smaller for the remaining 
values of r + . The values of £ considered lie above the critical values. 

Table H] contains the comparison of analytic versus numerical results for the two lowest 
QNMs at d — 6. The agreement is fairly good. Let us comment that the two critical values 
for £ read: £o = 0.90, £i = 1.86 for r + = 0.99. They are smaller for the remaining values of 
r + . The values of £ considered lie above the critical values. 

In figure [15] we show exact values of the imaginary part of the lowest purely dissipative 
mode for £ = 1 and £ = 2 versus r + for d = 5, along with a fit encoding the divergence 
at the critical point r + = 1. These results are successfully reproduced numerically as can 
easily be seen in figure [TBI The corresponding results for d = 6, where no analytical results 
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Table 3: Comparison of analytic vs numerical results for the two lowest QNMs at d — 5 
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Table 4: Comparison of analytic vs numerical results for the two lowest QNMs at d = 6 



are available, are contained in figure [TTj 
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Figure 15: Exact results for the imaginary part of the lowest purely dissipative mode for 
£ = 1 and £ = 2 v s r + for d = 5. 



5 Conclusions 



We have studied the perturbative behaviour of the topological-AdS black holes. We 
have calculated both analytically and numerically the QNMs of scalar perturbations of 
these black holes. 

Analytical calculations show that for small black holes at any dimension there is a 
critical point (at r + = 1) below which the real part decreases with n, having a positive 







-20- 



0.9 0.91 0.9 



0.94 0.95 0.96 0.97 0.9. 



Figure 16: Numerical results for the results for the imaginary part of the lowest purely 
dissipative QNM in five dimensions vs r + for £ = 1.0 (crosses) and £ = 2.0. 
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Figure 17: Imaginary part of the lowest purely dissipative QNM in six dimensions vs r + 
for £ = 1.0 (crosses) and £ = 2.0. 



slope, whereas above the critical point the oscillatory modes increase with a negative slope. 
We also found that below the critical point there is a critical value of £ below which there 
are purely decaying modes while above the critical point there are only oscillatory modes for 
any £. In five dimensions the QNMs of scalar perturbations of TBH-AdS can be obtained 
explicitly from the Heun function which solves the five-dimensional wave equation. 

These results are also supported by numerical investigations of the QNMs. The nu- 
merical results show clearly a change of slope of QNMs around a critical temperature for 
dimensions d = 4,5, 6. For larger dimensions the root finding algorithm is difficult to be 
implemented. This is connected with the observation that for higher dimensional theories 
the regularization needed for the convergence of the series becomes less efficient. 
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